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Abstract 



A tight binding representation of the kicked Harper model is used to obtain 

an integrable semiclassical Hamiltonian consisting of degenerate "quantized" 
^^ I orbits. New orbits appear when renorniahzed Harper parameters cross integer 

multiples of tt/2. Commensurability relations between the orbit frequencies 
^P ' are shown to correlate with the emergence of accelerator modes in the classical 

phase space of the original kicked problem. The signature of this resonant 
On . transport is seen in both classical and quantum behavior. An important 

feature of our analysis is the emergence of a natural scaling relating classical 



and quantum couplings which is necessary for establishing correspondence. 
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The search for evidence of complexity in quantum systems when the corresponding clas- 
sical system exhibits chaotic dynamics is an active frontier in nonlinear science. Kicked 
systems with classically chaotic dynamics have provided useful theoretical |]I],0] and exper- 
imental paradigms for elucidating the fundamental issues of quantum chaos. Systems 
where quantum states can be either extended or localized provide the added challenge of 
finding signatures of the localization-delocalization transition in the corresponding classical 
system. In this letter, we consider the kicked Harper system to establish precisely such a 
correspondence, within the context of fascinating issues such as anomalous transport Q and 
dynamical localization [Q]. 

The traditional approach starts with the classical analysis and then searches for signa- 
tures in the quantum behavior. By contrast, we begin with the quantum model and, using 
a novel semiclassical approach, predict a series of enhancements to transport in the system. 
These 'phase transitions' are shown to correspond to ballistic to diffusive transitions in both 
classical and quantum dynamics, and are due to classical accelerator modes. These modes, 
in turn, are known to result in anomalous transport [Q. Our methodology is applicable to 
a variety of kicked systems. However, we will illustrate our method in the kicked Harper 
model, 

H = Lcos{p) + Kcos{q)Y.6{t-k) , (1) 

k 

where the relative values of the parameters K and L determine the dynamics. The classical 
dynamics reduces to the two-dimensional area preserving Harper map, 

Pk+i = Pk + K sm{qk) 

Qk+i = qk - -Lsin(pfc+i) • (2) 

This mapping is analogous to the standard map or kicked rotor, with the dynamics becoming 
globally chaotic when both K and L cross critical values. However, neither system is ergodic 
as islands, arising via bifurcations, exist in the chaotic sea leading to the possibility of 
anomalous transport. In particular, hierarchical island chains for K = L m the vicinity of 
27r have been shown to result in fractional kinetics |Q. 

By contrast, the quantum dynamics of the kicked Harper model is quite different from 
that of the quantum kicked rotor. Most significant is the existence of both extended and 
localized states depending on K and L. Earlier studies indicate that extended states exist 
for K > L while localization occurs for i^ < L ||^, even though the K = L critical line, 
demarcating the onset of localization, has been rigorously established only in the absence 
of kicking. The existence of this transition allows the possibility of both ballistic as well as 
diffusive transport in the system. 

Our analysis emphasizes the dynamics along the line K = L and, using a novel semiclas- 
sical approach, predicts a rich variety of behavior. The key is the mapping of the temporal 
dynamics of the quantum kicked model to a continuum time-independent integrable Hamil- 
tonian, within which the complexity of the original problem manifests itself in terms of 
degenerate classical orbits. These orbits are indexed by an integer ra, where the n^^ orbit 
appears when the renormalized coupling K = K/{2h) = L exceeds mT/2. The phase space of 
this semiclassical Hamiltonian resembles the unkicked Harper model except that only certain 
orbits are allowed. Therefore, in this semiclassical formulation, kicking the system selects a 



discrete number of classical trajectories of the unkicked Harper system. These additional or- 
bits are shown to enhance both classical and quantum transport in the kicked Harper model. 
Finally, classical and quantum results display good correspondence provided the classical 
parameter is set to be K. This implies an h scaling considerably different from that used 
earlier |0. 

We start from the association of the time-dependent quantum Hamiltonian with a tight- 
binding lattice problem [^. This is obtained by projecting individual quasienergy states, 
associated with the single-kick evolution operator, onto angular momentum eigenstates \m > 
with eigenvalues p = mh. The projection coefficients Um satisfy 

Y^ Jr{K) sin [- (uj - 2L cos(m/i)) - —]um+r = 0, (3) 

where u is the quasienergy corresponding to the eigenstate. Following Wilkinson 0, we 
consider the continuum limit of the lattice model by defining Xm = ^rth and ip{xm) = i^m- 
Thus, {ipm+r+ipm-r) is replaced by 2 cos{rnp)i/j{xm) as p is the generator of space translation. 
Taking u = 0, the tight-binding model can be written as Hip{x) = 0, where the Hamiltonian 
operator H is, 

H= sin(Icos(x)) [MK) + 2^(-l)V2,(ir) cos(2rp)" 

+ cos(Z cos(x)) [2 Y.{-lY'^J2r-i{K) cos((2r - l)p)] . (4) 

Note that the lattice problem Htp{x) = can be viewed as an eigenvalue problem restricted 
to a Hilbert subspace defined by eigenvalue E = 0. The terms in the square brackets H 
can be summed exactly to give cos {Kcos{p)) and sin {Kcos{p)) respectively. Therefore, the 
resulting semiclassical Hamiltonian simplifies to 

Hsci = sin[K cos (p) + Zcos(x)] . (5) 

The fact that the Hamiltonian dynamics is constrained to E = implies that the classical 
phase space trajectories satisfy 

K cos{p) + Lcos(x) = ±mi, (6) 

where n is an integer. This can be viewed as the "quantization condition" for the phase 
space of the unkicked Harper equation given by Kcos{p) + Lcos{x) = E. Thus, within 
this semiclassical formulation, the kicking effectively quantizes the classical orbits. Another 
interesting feature, which we exploit later, is the fact that h appears only in K and L and 
not as an independent parameter. Changing h simply renormalizes the couplings K and L. 
Note also that Hgd is symmetric in x and p. 

The phase portrait resulting from Hgd consists of degenerate orbits indexed by the integer 
n. For K < 7r/2 and L < 7r/2, there exists only a single orbit and additional degenerate 
branches appear as K and L exceed threshold values K^ = n7i/2. The quantized orbits are 
not necessarily closed, unless K = L. For the rest of the discussion, we will confine ourselves 
to this symmetric case. In the corresponding quantum problem, the condition K = L for 
closed orbits suggests a phase boundary for a localization-delocalization transition |Q. For 
K < 7r/2 this boundary is the separatrix labeled by n = 0. Therefore, K = 7r/2, when the 



next orbit appears, can be viewed as the semiclassical threshold for the onset of diffusion in 
the kicked Harper rotor |^. 

The frequencies of the degenerate orbits are obtained from the equations of motion 
resulting from Hs^i 

p = i?(-l)"sin(x) 

x = -i?(-l)"sin(p) . (7) 

For any K, the frequency of the n^^ orbit is given, in terms of the complete elliptic function 
/C, by 

/„ = ir/[4/CK)], (8) 



where the modulus a^ = J[l — {K!^/ KY]. The turning points of these quantized orbits (for 



n > 0) a; = 2 cos (JK^/K) imply that, aX K = K'^, the new orbits appear at the origin 
(similarly, for n < they occur at (7r,7r)). At threshold, the n^^ orbit frequency f* = n/4 
is commensurate with the kicking frequency which was taken to be unity (see (|I])). As ^ 
increases, the frequencies increase monotonically and are, in general, incommensurate except 
at special parameter values. For example, for K ^ IAtc , we have /1//2 = 4/3. /i and /2 
are commensurate with the kicking frequency at several values: /i = 1/3 at ^ ^ .86tt, and 
/i = 1/2 at K ^ 1.687r while /s = 2/3 at K Ri l.Tlvr. The significance of these 'critical' K 
values will be seen shortly. 

In general, the quantum transport is expected to be enhanced by tunneling between 
these classically distinct orbits [0. Whenever a threshold value K^ is crossed, a new branch 
appears which may enhance the tunneling probability. We further expect that the transport 
may also be influenced by the commensurability relations between the competing orbital 
and kicking frequencies. 

Motivated by the possibility of phase transitions leading to the enhancement of diffusion, 
we now investigate classical as well as quantum transport in the original kicked Harper 
model. Some of the threshold values predicted by our semiclassical analysis are curiously 
reminiscent of known special values K = L = 2TTm beyond which accelerator modes appear 
in the Harper map. These are known to result in either ballistic or anomalous transport. 
On identifying K with the classical kicking parameter, the K* values include the thresholds 
for accelerator modes. This scaling emerges naturally from our semiclassical analysis and 
motivates the investigation of classical transport in the neighborhood of i^ = K^ = mT/2. 

Interestingly, a number of secondary fixed points of the classical map (eq.(|^)) appear 
at K^. Of these, a large number are accelerator modes, including the low-period ones 
associated with fractional kinetics [^]. For example, at ^ = 2n7r, four period-1 fixed points 
at (±7r/2, ±7r/2) appear onp = q and p = —q symmetry lines. Increasing K leads to pairs of 
period-1 points straddling the symmetry lines and, finally, hierarchies of island chains which 
result in anomalous transport. The parametric window for this behavior extends till K k, 6.6. 
For odd-multiples of tt {K = {2n + l)7c), analogous behavior with period-2 accelerator modes, 
one of which is (0, 7r/2) —>■ {tt/2,0) —>■ (0, 7r/2 + 2tiIk)i is seen. The island hierarchies at 
larger K are exemplified in Figure. |I|(b). Topological structures consisting of boundary layer 
island chains and, closer to the thresholds, the more exotic blinking and shrinking/blowing 
islands HS strongly suggest anomalous transport within parametric windows. Additional 



accelerator modes also appear at certain critical values, in particular, we have period-2 
modes at K ^ .STvr and K ?a l.SSvr on the p = q symmetry line. Note that these are close 
to parameters for the commensurability conditions /i = 1/3 and /1//2 = 4/3, respectively. 

The classical behavior at A' = (n + l/2)7r is more interesting with period-6 fixed points 
appearing on the symmetry lines p = 0,q = 0,p = q, and p = —q. As illustrated by the 
iterates {tt/2,tt/2) -^ (7r/2, {n + l)7r) -^ {{n + 1)tt, {2n + 3/2)7r) -^ {{2n + 3/2)7r, {2n + 
3/2)77) -^ {{2n + 3/2)7r, (n + l)7r) -^ ((n + l)7r,7r/2) -^ (7r/2,7r/2), these fixed points are 
non-accelerating despite leaving the fundamental interval [0, 27r] at intermediate steps. As 
K increases, these orbits also move off the symmetry lines as stable islands, as seen from 
Fig. |l|(a). Even though the fixed point is non-accelerating, an initial condition in the vicinity 
of the fixed point returns outside the fundamental domain after 6 iterations. In this sense, 
accelerated transport can occur at these K. 

The effect of these accelerator modes on classical transport is illustrated in Fig. |^(a) 
in terms of the exponent /i defined by < p^ >oc t'*, where the braces denote an ensemble 
average [1^. Many spikes in /i are seen as the parameter K is changed. The sharp, resonant. 



structures corresponding to ballistic transport, /i ^ 2, are due to the narrow ranges of 
stability of the accelerator modes. Interestingly, these peaks appear very close to the K* 
values as well as K values where the orbital frequencies exhibit commensurability. A few of 
these special values are indicated by the arrows in the figure. Thus our quantum analysis 
is providing information about the classical transport properties. Note that for smaller 
K/tt ^ 0.6, a stochastic web dominates the classical transport. 

As shown in earlier work @,|l3]! enhancement of the quantum localization length ^ can be 



used as a signature of classical ballistic behavior. Here we utilize this idea while, at the same 



time, contrasting the scaling suggested by our analysis with that used earlier p|,|6|,pJ]| . The 
resonant transport in classical and quantum behavior shows good correspondence provided 
the classical stochasticity parameter is taken to be ^ = K/{2h). This is quite different 
from Kg = Ksin{h/2)/[h/2) obtained earlier [^ by contrasting correlations in the classical 
and quantum dynamics. K emerges naturally from a 'thermodynamic' formalism, as the 



tight-binding (eigenvalue) problem deals with only a single quasienergy state |]T2|. As such, 
this may be valid beyond a certain time scale. We consider moderate values of h and, as 
seen from panels Fig. |](b) and(c), spikes in ^ appear close to classical values only when K, 
and not Kg, is considered. At these values of h, the issue of reaching large t behavior does 



not arise. Also, as in earlier work |TT[, the quantization scale is much larger than the size 
of the classical structures causing the accelerated transport. The intriguing possibility, for 
smaller h, of a cross-over from Kg to K at longer times warrants further exploration. 

Figure ^ shows the quantum lineshape at one of the peaks (a) contrasted with the shape 
away from the peak (b). It is clear that the quantum transport at the peak has a ballistic 
component which results in the well-defined shoulders in Fig. ^(a) as compared with panel 
(b). 

Despite the unphysical nature of the Harper system, resonant transport at K* and other 
critical values can be tested experimentally. As the origin of resonances has been related to 



accelerator modes, analogous behavior exists in the kicked rotor (standard map) 0,0 and 
oscillator (web map), both of which are experimentally realizable in atom optics. In fact, 
even the question of scaling as discussed above could also be resolved experimentally in this 
context. 
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FIGURES 

FIG. 1. Phase space portraits of the accelerator islands beyond the thresholds K = K^ = wk jl. 
(a) K = 5.04(n = 3) where one of the six islands is shown; (b) K = 3.199 (n = 4) where only one 
of four symmetric structures is shown. Note the size of these structures in relation to the h values 
considered. 

FIG. 2. (a) Variation of the classical transport exponent /i with the classical stochasticity pa- 
rameter set to K. An ensemble of 1001 initial conditions along the line p = was iterated 10^ 
times, and the last 5 x 10*^ iterations were used to obtain /i. The arrows at the top indicate some of 
the special values arising from commensurability and thresholds K^ (see text), (b)-(c) Variation of 
the localization length ^ with K (solid) and Kg (dashed), h = 2iTa^ and (b) M = 1 and (c)M = 3, 
where a refers to the golden mean. A fixed interaction time of 4000 kicks was considered. 

FIG. 3. Quantum hneshape at the end of 4000 kicks. (a) K = 6.4402648 while (b) 
K = 6.3146012. % = 27Ta^ and case (a) is at a peak in Figure 2(c) while (b) is off the peak. 
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